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We consider the excitation dynamics of the two-photon 1S-2S transition in a beam of atomic 
hydrogen by 243 nm laser radiation. Specifically, we study the impact of ionization damping on 
the transition line shape, caused by the possibility of ionization of the 25* level by the same laser 
field. Using a Monte-Carlo simulation, we calculate the line shape of the IS -25 transition for 
the experimental geometry used in the two latest absolute frequency measurements (M. Niering 
et al, PRL 84, 5496 (2000) and M. Fischer et al, PRL 92, 230802 (2004)). The calculated line 
shift and line width are in excellent agreement with the experimentally observed values. From this 
comparison we can verify the values of the dynamic Stark shift coefficient for the 15 -25 transition 
for the first time on a level of 15%. We show that the ionization modifies the velocity distribution of 
the metastable atoms, the line shape of the 15 - 25 transition, and has an influence on the derivation 
of its absolute frequency. 

PACS numbers: 32.70.-n, 32.90.+a, 42.50.Hz 



I. INTRODUCTION 

High- precision spectroscopy of the 15-25 transition 
in atomic hydrogen and deuterium provides an essential 
contribution to the determination of the Rydberg con- 
stant, the Lamb shift [l|, the 25 hyperfine interval [2J as 
well as it allows to study nuclear properties of the proton 
and the deuteron Recently it has been shown, 
that monitoring the absolute frequencies of narrow clock 
transitions (like the 15-25 transition in hydrogen and 
quadrupole transitions in the Hg + and Yb + ions) over 
a prolonged time interval, opens an opportunity to set 
a stringent restriction on the possible slow variation of 
the fine structure constant |5|, |6J , testing very fundamen- 
tal aspects of physics. This laboratory approach relies 
on the accuracy of the 15-25 frequency measurement, 
which has now reached the level of 1.4 x 10~ 14 , but in 
principle has still room for improvement, since the Q- 
factor of this transition is about 2 x 10 15 . The accuracy 
which is achievable experimentally is restricted by a num- 
ber of systematic effects which shift the frequency and 
change the line shape. So a further increase of accuracy 
requires not only an upgrade of the experimental setup 
and an improvement of the signal-to-noise ratio, but also 
a tighter control of the systematic effects. Among these, 
the dynamic Stark shift, the second order Doppler ef- 
fect and the time-of-flight (TOF) broadening have been 
considered as the most important issues until recently. 

In this article, we analyze how the presence of ion- 
ization losses of the 25 state via the absorption of 243 
nm photons contributes to the line shape of the 15-25 
transition in the experimental configuration used for the 
most recent absolute frequency measurement 0,0- Com- 



paring the experimentally obtained intensity-dependent 
line width and frequency shift of the 15-25 transition 
with results of a numerical simulation, we can also in- 
directly validate a number of atomic parameters like the 
two-photon transition matrix element, the AC Stark shift 
coefficient and the ionization coefficient, which were re- 
cently re-derived in a unified formalism H|. In addition, 
Monte Carlo simulations allow to study systematic effects 
which are difficult to quantify experimentally. From this 
we gain new insights on the origin of the scatter in the 
15-25 absolute frequency data observed in the experi- 
ments of 1999 and 2003 [EGi 

The article is organized as follows: In Sec. II, we de- 
scribe in brief the experimental apparatus emphasizing 
points important for the current analysis, while the nu- 
merical model is given in Sec. III. We analyze the results 
in Sec. IV, followed by the conclusions. 



II. EXPERIMENTAL SETUP 

The basic experimental setup used for the spectroscopy 
of the 15-25 transition in a cold beam of atomic hy- 
drogen is described in detail in Q. In this section we 
will focus on a description of the apparatus depicting the 
main features necessary for the current analysis and will 
present the relevant experimental results. 

The dipole- forbidden 15-25 transition is excited by 
means of a Doppler-free two-photon driving scheme in 
counter-propagating waves of equal frequencies. The nec- 
essary 243 nm radiation is produced by doubling the cw 
light emitted by a 486 nm dye laser in a /3-barium borate 
crystal. The laser is stabilized to an external ultra-stable 
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FIG. 1: Excitation of cold hydrogen atoms by 243 nm light in 
the enhancement cavity. Atoms escaping the cold nozzle are 
collimated by two diaphragms Dl and D2 before entering the 
detection zone. All parts except the chopper assembly and 
the uv photodiode are placed in vacuum of l(T 7 -lCr 6 mbar. 
The diameters of the nozzle and the diaphragms are drawn 
not to scale and are much smaller as suggested by this graph. 

cavity and has a spectral line width of about 60 Hz while 
the drift of the laser frequency is less than 0.5 Hz/s 
The frequency of the 486 nm laser is measured using the 
frequency comb technique 0, which is described in 
detail in |nj]. 

To efficiently excite the 1S—2S transition in atomic 
hydrogen we enhance the 243 nm radiation in a linear 
cavity consisting of a plane incoupling mirror (transmis- 
sion Ti = 2.4 x 10~ 2 ) and a concave outcoupling mirror 
(radius R = — 4 m, T 2 — 1.36 x 10~ 4 ) separated by a 
distance of 29 cm as depicted in Fig. Q The cavity loss is 
dominated by the diaphragms Dl and D2 that partially 
touch the laser mode so that the input coupling mirror 
roughly corresponds to impendence matching condition. 
The radiation leaking out of the cavity is monitored by 
a silicon photodiode which is periodically calibrated by a 
power meter (Fieldmaster, Coherent Inc.). The photodi- 
ode readout has been averaged over one second which is 
the laser dwelling time when recording the spectra. Di- 
viding the photodiode readout by the transmission of the 
outcoupler T 2 we get the total power circulating in the 
enhancement cavity per direction. The averaged power 
is recorded for each measured spectrum. 

Atoms produced in the microwave gas discharge escape 
a cold nozzle (T ~ 5 K) of 1.2 mm in diameter and enter 
the interaction zone of length £' = 15.0 cm, restricted 
by two diaphragms Dl and D2 of respective diameters 
2rDi = 1.3 mm and 2r^ 2 = 1-4 mm separated by £ — 13.6 
cm. During the flight through the Gaussian profile of the 
cavity TEMoo mode, a part of the atomic hydrogen is 
excited to the metastable 2S state. Excited atoms are 
detected by quenching in a dc electric field through the 
release of a Lyman-alpha photon @ . 

To set an upper limit for the velocity of the atoms 
that contribute to the signal and thus to reduce the 
influence of the velocity-dependent systematic effects, 
a time-of-flight detection technique [3, M, is used. All 
atoms with velocities higher than i> max escape the detec- 
tion zone before the start of the detection and the con- 
tribution of these atoms is excluded from the resulting 
time-delayed spectrum. The cut-off velocity is defined as 
v max = V '/At, where At is a delay between closing of 
the 243 nm laser light and the start of detection. 



The 486 nm laser is scanned stepwise across the reso- 
nance and for each laser frequency a multichannel scaler 
records the number of detected Lyman-a photons for 12 
different delays An = 10, 210, 410, . . . , 2210 fis in the 
corresponding time windows {At, 3 ms}. 

Using this detection technique we simultaneously 
record 12 lines, each containing the contributions of 
atoms in the respective velocity range < v < v max {Ari) 
(i = 1, ...,12), selected from the same initial velocity 
distribution f(v). As follows from the velocity dis- 
tribution in a one-dimensional thermal atomic beam ef- 
fusing from thermalizcd gas volume through a hole in a 
thin wall is given by: 

f(v)^(v/v ) 3 eM~(v/v ) 2 }, (1) 

where vq — ^j2ksT /toh is the most probable thermal 
velocity for a given temperature T and mn is the mass of 
hydrogen atom. As discussed in Ref. 0] and references 
therein, expression (JJJ is valid for perfectly collimated 
beams and should be further modified in other cases also 
because the thin wall assumption is not fulfilled in our 
case. The collimation angle of the atomic beam in our 
experiment is about 0.01 which justifies the use of the 
mentioned approximation. The impact of this approxima- 
tion on the 1S-2S frequency shift will be discussed in 
Section ITVCl 

The increase of the delay At reduces the maximum 
velocity !i max (Ar) of atoms contributing to the signal 
and decreases the Lyman-a count rate. Correspond- 
ingly, the contributions of the time-of-flight broadening 
and the second-order Doppler effect are reduced, which 
in turn leads to a narrowing and a symmetrization of 
the 1S-2S transition line shape. Hence, at the cost of 
a lower signal-to-noise ratio, lines with lower systematics 
are recorded. Using this experimental approach and the 
line-shape model |!j one can correct for the second order 
Doppler effect which is necessary for the accurate deter- 
mination of the absolute 1S — 2S transition frequency. 

In this paper we will study three types of intensity 
related line broadening effects: Because the 2S state is 
close enough to the continuum that a single additional 
243 nm photon is sufficient to ionize the hydrogen atom 
reducing the effective interaction time. In addition the 
spatially dependent AC Stark shift and the usual power 
broadening caused by Rabi flopping are present. 

For simplicity, we will analyze only hydrogen lines 
recorded at a particular delay, At = 1210 fis which corre- 
sponds to f max — 120 m/s: this is a compromise between 
simple line shape and sufficient signal-to-noise ratio. At 
this delay the line shape can be well approximated by a 
Lorentzian and the maximum second-order Doppler ef- 
fect does not exceed 200 Hz in modulus, the typical line 
width being 1-2 kHz at 121 nm. Therefore we consider 
contributions to the width of a Lorentzian line, which 
provides a much clearer insight into the nature of the dif- 
ferent broadenings than the complete line shape model 
for all At. 
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power per direction P [W] 



in the following that this coefficient is formed by three 
important contributions, namely power broadening, in- 
homogeneous AC Stark shift in a Gaussian beam profile 
and the influence of the ionization of the 25 level. 

In the next section we will analyze the structure of 
these coefficients by the help of a numerical analysis per- 
formed for an ensemble of hydrogen atoms excited in the 
given experimental geometry. 



III. EXCITATION OF A BEAM OF ATOMIC 
HYDROGEN IN THE PRESENCE OF 
IONIZATION 
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FIG. 2: Results from 59 line scans: (a) center frequencies 
and (b) line widths (FWHM) both derived from Lorentzian 
fits versus the excitation power circulating in the cavity per 
direction, P. The power of 1 W per direction corresponds to 
a total on-axis focal intensity of 16 MW/m 2 averaged over 
one wavelength of the standing wave. 



To correct for another important systematic effect, the 
dynamic Stark shift, we recorded a number of 15-25 
spectra at different excitation powers during each day of 
measurement Here we use the same set of data for 
the analysis. 

All lines recorded at At = 1210 /is have been fitted 
by Lorentzian functions and then the absolute frequen- 
cies of the line centers v\s-2S and line widths r to t were 
plotted versus the excitation light power P (per direc- 
tion) . The results for the data set recorded on one of the 
days of measurement and the linear extrapolation to zero 
intensity are presented in Fig. [5] as an example. 

The intensity-dependent line shift v\s—2s{P) is mainly 
caused by the dynamic Stark shift |8( with small contri- 
butions of velocity-selective ionization. As in Refs. 0,0 
we fit the data with a linear function (since for low in- 
tensities the AC Stark shift for 15-25 transition is pro- 
portional to the intensity), but instead of extracting just 
the absolute frequency at zero power we also determine 
the slope of the fit, fc s hiftj as shown in the figure. Corre- 
spondingly, we consider the slope of the linear fit of the 
rtot(-P) dependence and call it Abroad- It will be shown 



The calculation of the 25 excitation dynamics is based 
on the equations of motion for the matrix elements 
of the density operator p' in the rotating wave ap- 
proximation |8|]. The laser field of angular frequency 
uil ~ 2ir x 1233 THz, corresponding to A = 243 nm, and 
intensity I excites the two-photon transition between the 
ground state \g) (15) and the excited state |e) (25). The 
frequency detuning is defined to be the difference between 
the doubled laser frequency and the resonance frequency 
Auj = 2lul — uj' eg , where hoj' eg is the energy difference 
between the 15 and 25 states of a perturbed atom (see 
Eq. iJSJ), ignoring the the hyperfine structure 0. 

In the experiment, the excitation field intensity I de- 
pends on the position inside the excitation region, be- 
cause the laser profile is of Gaussian type. The fact that 
each atom flies along a certain trajectory with an indi- 
vidual speed through this intensity profile, renders the 
intensity a time dependent quantity I(t) in the frame of 
the travelling atom. Because the dynamic Stark effect is 
proportional to the intensity, also the total detuning as 
measured by the atom, is a function of time, Aw(t). In 
addition the upper state |e) decays to the continuum due 
to photoionization with a time-dependent rate propor- 
tional to the laser intensity. One can write the equations 
of motion for the density matrix, neglecting spontaneous 
decay of the 25 state [8j . After factoring out the fast op- 
tical oscillations from the components of the density ma- 
trix p via p' ~ pe~xp(—2iuji,t) we obtain the two-photon 
Bloch equations: 



= -n(t)lm(p'(t)), 



(2a) 



= -iAw(t)p'(t) + i-¥-(p' (t) - pfjfi) 



d_ 

dt 



-^PU*), (2b) 
p' ee (t) = n(t)Im(p'(t))-7i(t)p' e M (2c) 



where the time-dependent two-photon Rabi frequency 
is fi(i) = 2(27r/3 ge )(m e //z) 3 7(t). Here m e is the elec- 
tron mass, p the reduced mass, and the value of j3 ge 
is taken from (3 ge = 3.68111 x 1(T 5 Hz(W/m 2 )- 1 
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FIG. 3: (a) - normalized responses of a single atom illuminated by a 1 ms rectangular pulse with the intensity of I = 4 MW/m 2 
with ionization channel (71 7^ 0) and without it (71 = 0). For simplicity, neither the AC Stark shift nor the Doppler effect 
are taken into account, (b), (c) - excitation probabilities for atoms flying along different random trajectories through the 
Gaussian intensity profile corresponding to P = 300 mW. (b) - calculations for 50 randomly chosen trajectories in presence of 
the ionization channel, (c) - calculation for the same ensemble of trajectories without ionization taken into account. 



for the 1S-2S transition. In the same paper the ioniza- 
tion rate is given by Ji(t) — 27r/3i on i(e)(m e //Lt) 3 /(i), where 
the ionization coefficient reads /3i on i(e) = 1.20208 x I0~ 4 
Hz (W/m 2 )- 1 . 

Experimentally an ensemble of atoms contributes to 
the signal, with each individual atom having a different 
speed and trajectory through the interaction zone. We 
take into account only those trajectories, which are con- 
fined by the diaphragms Dl and D2 (see Fig.Q. It is as- 
sumed, that the starting points on the diaphragm Dl and 
the exit points on the diaphragm D2 defining each tra- 
jectory are randomly and independently distributed over 
the corresponding diaphragm with constant area density. 
It means, that the angle tp > between the cavity axis 
z and each individual trajectory is also randomly and 
uniformly distributed for r^, rrj2 <C I'- the atom can in- 
tersect the excitation light beam as well as fly parallel to 
its axis. The velocities are also chosen at random accord- 
ing to a Maxwellian distribution like in Eq. at T = 5 
K, as in the experiment the day-averaged readouts of the 
temperature sensor, placed on the cold nozzle, ranged 
between 4.95 K and 5.25 K. 

After this longitudinal velocity seeding, all atoms 
which would reach the detector earlier than the time 
at which the delayed detection started (that is before 
At = 1210 /is), were rejected from the ensemble. For 
each remaining atom, one can then easily calculate the 
time-dependent intensity I(t), as well as the total in- 
teraction time. The laser is modelled by a Gaussian 
beam profile with parameters defined by the experimen- 
tal geometry. The radial intensity distribution is there- 
fore given by 7(r) = (AP/ttw 2 ) exp[— 2r 2 /w 2 ], where 
w 2 = Wq(1 + {zX/ttw 2 ) 2 ). The beam slightly diverges 
along the z-axis, and the beam waist on the incoupler 
plane equals wq = 283 /im. Here P is the power per 
direction, as used in Fig. [5] 

It is necessary to point out, that atoms flying through 
a standing wave in the cavity interact with a highly in- 



homogeneous laser field consisting of nodes and maxima 
with the periodicity of 121.6 nm in z-direction. But even 
for a very slow atom travelling at a longitudinal veloc- 
ity of v = wo/300 ~ 1 m/s, Eq. QJ, the intensity will 
be modulated at a frequency of about 10 GHz, which is 
6 orders of magnitude higher than the typical Rabi fre- 
quency f2, which governs the population dynamics. Thus 
for the calculation of the density matrix elements, we can 
equivalently use a laser field with an averaged intensity, 
in a similar way as in applying the rotating wave approx- 
imation. 

The general detuning Aw(t) along the trajectory can 
now be expressed as 

Aw(i) = 2u L -u' eg = (3) 

= 2u L - LO eg + (-) -27T/3 ac /(<), 

where hio eg is the energy between the unperturbed 15 
and 25* states. The other summands take into ac- 
count the second order Doppler effect and the intensity- 
dependent frequency shift caused by the real part of the 
dynamic Stark effect. In turn, the imaginary part of the 
dynamic Stark effect is related to ionization. According 
to@,/? ac = [p ac (2S)-(3 ac {lS)](m e /fi) 3 = 1.66982xl0- 4 
Hz (W/m 2 )- 1 . 

The time evolution of Eqs. J3J) for a specific atom and 
its corresponding interaction time is calculated numeri- 
cally for each given laser frequency u>l which is scanned 
stepwise across the resonance. The resulting value of p' ee 
is stored together with wt to define the contribution of 
each individual atom to the line profile. 

Fig-El 1 illustrates the spectral response of a single atom 
at rest to the excitation with the rectangular pulse of 
243 nm radiation. This case is analyzed in detail in 
One observes, that for moderate excitation intensities of 
a few MW/m 2 and for short interaction times less then 
1 ms, the contribution of the ionization broadening to 
the total line width is practically negligible. Indeed, for 
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a 1 ms rectangular pulse of intensity 1 = 4 MW/m 2 the 
ionization broadening is a very small effect, though the 
ionization rate itself in this case equals 7; = 480 Hz. 

Fig-El 3 illustrates how individual atoms flying through 
the excitation region contribute to the observed line 
shape. Depending on the individual trajectory and the 
velocity of an atom it will be exposed to a time vary- 
ing intensity. The figure shows the resulting individual 
line shapes, over which the detector integrates to yield 
the observed line profile. Note the variations in ampli- 
tude, shift and width of the individual lines due to vary- 
ing speeds (second order Doppler shift), interaction times 
(excitation and ionization probability, TOF broadening) 
and sampled intensities (excitation and ionization prob- 
ability, AC Stark shift, power broadening). 

In contrast to the physical case given in Fig. EJd with 
all effects taken into account, Fig. EJ; shows the excita- 
tion probabilities for the same ensemble of 50 trajectories 
calculated without ionization. Apparently the ionization 
preferably removes the slow atoms from the initial dis- 
tribution, because they have more time to interact with 
the field. But these atoms would have contributed the 
strongest and the narrowest spectral lines to the resulting 
collective spectrum if no ionization had taken place. As 
a result, the detected line becomes weaker and broader. 

In the results section we will discuss the simulated line 
profiles in more detail and compare the results with ex- 
perimental data. 

To perform the described Monte-Carlo simulation we 
used the package "Mathematica 5" 0] . The simulation 
has shown good convergence and an arbitrary ensemble 
of 10 000 atoms has been used to calculate each transition 
line shape. The input parameters used by the program 
can be divided into four groups: (i) atomic structure, (ii) 
geometric parameters (wq, tdi, td2, £, £' , nozzle radius), 
(iii) experimental conditions (P, T, At) and (iv) two 
program switches which allowed to artificially set 7; = 
(no ionization) or /3 ac = (no AC Stark shift). The last 
option allows to analyze separately the contributions of 
ionization, AC stark shift and power broadening to the 
line shape. 

Including all contributions to p' ee {Auji,), we obtain the 
simulated line shape, which we fit with a Lorentzian pro- 
file exactly like in the evaluation of the experimental data 
for the purpose of this study. In our simulations, we var- 
ied the excitation power up to P = 1.2 W per direction 
which covers and exceeds the whole experimental range. 



IV. RESULTS 

In this section the results of the numerical simulation 
will be compared with experimental results obtained in 
the measurement . We will analyze the frequency shift 
and the line shape of the 1S-2S transition and discuss 
possible origins of the excess data scattering observed in 
the frequency measurement experiments. 

We would like to point out here, that both the fre- 



quency shift and the line broadening are very small ef- 
fects of the order of 10 -12 compared to the optical carrier 
frequency. Only the improvements of the laser stability 
and the narrowing of its spectral line width along with an 
improved absolute frequency measurement have made it 
possible for the first time to perform a quantitative anal- 
ysis of the line shift and especially of the line broadening 
as a function of excitation intensity. 

Figure ^ illustrates the results of the different simu- 
lations (a, b) and their comparison with experimental 
observations (c, d). To evaluate the uncertainty of the 
simulation for 10 000 atoms due to the random elements 
of the simulation, we performed the simulation for 3 dif- 
ferent arbitrary ensembles of atoms for the same input 
program parameters. The variation is on the level of few 
per cent which is sufficient for this analysis. 

A. Line Center Shift 

As follows from Fig. the main origin of the fre- 
quency shift is the AC Stark shift. The physical case is 
depicted in a solid curve (71 7^ 0, /3 ac 7^ 0). The line shift 
reveals an observable non-linearity at higher powers. It 
could be explained by the following contributions: (i) the 
inhomogeneous AC Stark shift, (ii) the influence of ion- 
ization and (iii) the saturation of the transition. In the 
current analysis we took the non-linearity into account by 
fitting the simulations in Fig. 0Ji with a parabolic func- 
tion. 

For comparison with the full treatment, we performed 
the same set of calculations for three non-physical cases: 
with the ionization channel switched off (71 = 0, /3 ac ^ 
0), the AC Stark shift is switched off (7 ^ 0, /3 ac = 0), 
and both these effects are switched off (7; = 0, (3 ac = 0). 
We see, that the main contribution to the frequency shift 
is due to the real part of the dynamic Stark shift. The 
ionization, which reduces the number of atoms undergo- 
ing the highest perturbation in the intensive laser field. 
It also influences the frequency shift, but only if the dy- 
namic Stark shift is present as well. All data sets in 
Fig. EJl converge at zero power to the same value of 
—20(1) Hz which corresponds to the mean second order 
Doppler shift for the given beam temperature of T = 5 K 
and the delay time At = 1210 fj,s. This value is consid- 
erably smaller than the maximal shift of —200 Hz calcu- 
lated for u max (Ar). This difference is caused by the fact 
that the velocity distribution of atoms which contribute 
to the signal f'(v) is different from the velocity distri- 
bution of atoms in the beam f(v): neglecting additional 
losses it is given by 

f'(v) oc f(v)(l-v/v majc (AT)), v<v max , (4) 
f'{v) =0, v > v max . 

Consequently, the Doppler correction must take this 
weighting into account. A further correction to the ef- 
fective second order Doppler effect is invoked by the pre- 
ferred removal of slow atoms. As seen in Fig0^ (for 
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FIG. 4: Simulated dependencies of the frequency shift (a) and the line width (b) of the 1S-2S transition vs. excitation power 
per direction, P and a detection delayed by At = 1210 fis. The set of data marked /3 ac 7^ 0, 71 7^ corresponds to the 
physical case. The other simulations clarify the separate contributions of ionization, inhomogeneous AC Stark shift and power 
broadening (e.g. the data set with 71 = 0, /3ac = in (b) corresponds to pure power broadening). Note that the contributions 
cannot be added in any simple algebraic way. The data sets are fitted by parabolic functions. On the left side of (c) and (d) 
the the day-averaged values of fcshift and fcbroad with corresponding uncertainties are presented. The total 12-day average of the 
data is indicated by the dashed lines and the gray bars for the uncertainties. Simulations for two different beam temperatures 
(5 K and 4 K) and for different radii of the entrance diaphragm Dl (roi, t"di/2 and wo) are presented on the right part of the 
lower plots. 



ji 7^ 0, /3ac = 0) this preference and thus the associated 
shift is vastly independent of laser power for our experi- 
mental conditions. As it is impossible in our experiment 
to measure the Stark coefficients /3 ac and /3; on i separately, 
we analyze the values which can be derived both from 
the experiment and our numerical analysis, namely the 
slopes fcshift and fcbroad of the power dependencies given 
in Fig. and Figs. EJi, b. 

On the experimental side, we derive these values for 
each of the 12 days of measurement exactly in the same 
way as depicted in Fig. |2 The results for each day are 
presented in Figs. 0]c,d. One observes a significant scat- 
ter of the day-to-day data, with Xr — 3 (the reduced \ 2 
or Birge ratio |14|). which could be partly ascribed to er- 
rors on the power axis, coming from both the calibration 
of the photodiode measuring the 243 nm power and from 
temporal variations of the transmittance of the curved 
mirror in Fig. ^ Therefore we average these data with- 
out taking into account the individual uncertainties and 
obtain fc s mft = 1-61(6) Hz/mW and fcbroad = 2.25(11) 
Hz/mW. 

Extracting the coefficient fc s hift from the results of the 



numerical analysis, depicted in Fig. is not straight- 
forward, since the curve for the shift of the line center 
are not linear and therefore fcshift is not strictly defined. 
However, since the scatter in the experimental data is 
larger than the calculated nonlinearities, it is enough on 
the present level of accuracy to replace the parabolic fit 
by a linear fit in the case of the frequency shift-intensity 
dependence (fcshift), in order to evaluate the results in 
the same way as in the experiment. To reduce the error 
introduced by this procedure, we restrict the data points 
taken into account by the fit to those which fall into the 
experimental range of powers: 0.1 W < P < 0.5 W. The 
value resulting from the Monte Carlo analysis is denoted 
as fc*hrft r i n analogy to the notation in Fig. ^ 

To test the sensitivity of the model to the experimen- 
tal conditions and to evaluate the corresponding uncer- 
tainties we also performed a set of simulations for differ- 
ent beam temperatures, varying between 4 K and 6 K, 
slightly different interaction lengths and a reduced ra- 
dius n of the diaphragm Dl. All these tests show a low 
sensitivity of results to these parameters which allows 
to compare them with the experimental data. We have 
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performed calculations of fc^ST f° r different realistic ex- 
perimental parameters. The results of four selected cases 
are given in Fig. with an atomic beam temperature 
of 4 and 5 Kelvin and for a diaphragm size matched and 
centered to the laser beam to simulate the nozzle covered 
with a molecular hydrogen film. 

The last case corresponds to a possible situation, where 
a film of molecular hydrogen symmetrically covers the 
nozzle from inside in such a way that the nozzle radius 
equals either tdi/2 or wq = 283 /xm (the film eventu- 
ally touches the laser beam). The frozen nozzle radially 
restricts the atomic trajectories on the entrance to the 
interaction zone, which results in changing the intensity 
profiles. The growth of the H2 film was a characteris- 
tic feature of the hydrogen spectrometer working at the 
lowest beam temperature. The measurement cycle lasted 
about 15 min, after which the nozzle was heated to 20 K 
to melt the film. 

Comparing the averaged value fc s hift with the results of 
the numerical simulation, we can establish an agreement 
on the level of about 15% uncertainty with respect to the 
absolute values. As follows from the previous analysis, 
the intensity-dependent frequency shift is mainly caused 
by the AC Stark shift which contributes to the general 
detuning © ■ Using the linearity of the effect at low pow- 
ers one can conclude, that this analysis allows to verify 
/? ac at the same confidence level of 15% and therefore we 
can experimentally verify the results of the calculations 
of &.c = /? ac (2S) ~ Presented in @. 



B. Line Width 

The analysis of the line width is more complicated 
in comparison with the analysis of the frequency shift 
since a number of different effects contribute signifi- 
cantly to the line broadening. The line width consists 
of intensity- independent and intensity-dependent parts. 
The intensity-independent part includes the laser spec- 
tral line width and the time-of-flight broadening, while 
such effects as power broadening, inhomogeneous AC 
Stark broadening and ionization broadening contribute 
to the latter part. In addition, temporal intensity fluctu- 
ations of the laser radiation in the excitation region also 
induce an additional intensity-dependent broadening due 
to the variable Stark shift. Thus, assuming for simplicity 
the Lorentzian profiles, one can write: 

r t ot(P) = r(p) + (5) 

where T tot (P ) is the total line width, r to f is the contribu- 
tion of the time-of-flight broadening, and T(P) includes 
all intensity-dependent effects with T(0) = 0. The laser 
line width ri ascr carries a factor of 4 because we mea- 
sure it at 486 nm whereas all other transition frequencies 
correspond to 121 nm. It has been shown experimen- 
tally, that the spectral line width of our laser system 



doubles for the frequency doubling process. We under- 
stand that each of the contributions listed above has a 
complex spectral shape which can significantly differ from 
Lorentzian and the expression (JjJJ can be used only for 
rough estimations. Still, the simulations show that in our 
case the time-of-flight broadening has near-to-Lorentzian 
spectral line shape while heterodyne experiments indi- 
cate the same for the spectral line of our laser system. 

As is evident from Fig. 0Jd, the intensity-dependent 
part of the line width is formed by three different pro- 
cesses: (i) the power broadening, (ii) the inhomogeneous 
AC Stark shift in the laser beam profile and (iii) the 
ionization of the 2S state. The pure power broadening 
corresponds to the case (7; = 0, (3 ac = 0) and is due to 
Rabi flopping between the ground and the excited states, 
effectively reducing the lifetime of the population in the 
excited 2S state. The inhomogeneous AC Stark shift in 
the case (71 = 0, /3 ac 7^ 0) originates from a Gaussian 
intensity profile in the excitation region. The last signifi- 
cant broadening comes from ionization which contributes 
in a highly nontrivial way illustrated in Figs. |3J),c. 

The excessive loss of slow atoms has been taken into 
account in the line form model 0, Q as an empirical 
modification of the low-velocity wing of the Maxwellian 
w 3 -dependency in Eq. ■ There it was assumed that the 
slow atoms are missing from the beam because of colli- 
sions with the fast atoms (Zacharias effect). However, 
now we conclude that no such mechanism is necessary 
to fit our data and the selective ionization is the most 
significant effect responsible for this loss. Strictly speak- 
ing, one deals with three distinct velocity distributions: 
(i) the atoms emerging from the nozzle, (ii) the detected 
2S atoms, and (iii) the ground state atoms. Whereas the 
first distribution is assumed to be Maxwellian (Eq. Q), 
the second is relevant for determining the line shape. 

To compare with the experimental results, we also eval- 
uated the contribution of the temporal intensity fluctu- 
ations. In principle, this effect is equivalent to the con- 
sidered inhomogeneous AC Stark shift, but integrated 
over time, not space. In the experiment, we observed 
fast intensity fluctuations on the level of ±5% during the 
excitation cycle. We embedded this process into the sim- 
ulation by adding white noise with a 5% band to I(t). 
The effect lies within the uncertainty of the calculations 
for the given number of atoms used in the simulation. 

For the line width as a function of laser power, the 
slope is derived from the fit similar to the case of Fig.QJi. 
The value resulting from the Monte Carlo analysis is de- 
noted as Abroad- The line width, extrapolated to zero 
power, is r^f° r = 550(5) Hz which corresponds to the 
theoretical value for the averaged time-of-flight broaden- 
ing for the given delay r = 1210 /is in the absence of 
other broadening contributions (Fig^J)). 

The value fcbroad can be tested in the same way as 
fcshift with an uncertainty of 18% (see Fig. Due to 
the complex structure of this value and the number of 
contributions of different nature which add up to the to- 
tal line broadening, it is impossible to disentangle them 
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and to set separate restrictions on the transition matrix 
element and the ionization cross-section. Nevertheless, 
if we switch off any of the important broadening mech- 
anisms, as illustrated in Fig. the simulation will no 
longer agree so well with the experimental observation. 
This means that both the calculations performed in |£| as 
well as the current analysis adequately describes the ex- 
citation dynamics of the hydrogen atom in the presence 
of photoionization. 

For consistency of our analysis, we also performed an 
evaluation of the laser line width ri ascr . Studies of the 
spectral laser line width at 486 nm, based on the investi- 
gation of beat signals, resulted in a value of 60 Hz ^] . In 
the framework of this paper, we can make an estimate of 
its spectral line width with a different, indirect method. 
To this end, we compare the calculated line width at zero 
laser power T^f or with the experimentally observed line 
width r cxp , where the latter is obtained by averaging 
the extrapolated day values for the line width similar to 
Fig. resulting in r cxp = 775(20) Hz. From Eq. © 
we obtain ri aser = 56(5) Hz, which is consistent with the 
previous independent measurement |15| . The splitting of 
the (15, F = 1, m F = ±1) -> (25, F' = 1, m' F = m F ) 
magnetic components is negligible in the experiment. 

An attempt to directly evaluate the expected Lyman-a 
count rate in the given experimental geometry and thus 
test the theoretical values more directly fails due to the 
large uncertainties in the beam density, the fraction of 
atomic hydrogen in it and the detection efficiency. The 
presented indirect analysis happens to be much more sen- 
sitive to the absolute values of the transition matrix el- 
ements and the ionization cross section than the direct 
evaluations. 



C. Impact on the 1S — 2S Absolute Frequency 
Determination 

At the end of this section we discuss the conclusions 
of the current analysis regarding the absolute frequency 
measurement of the 1S — 2S transition in atomic hydro- 
gen. 

First of all, we want to point out the non-linearity of 
the frequency shift of the line center as a function of 
the laser power (Fig. ^i). Our analysis shows, that our 
previous linear fit leads to an error of 10-20 Hz in the 
determination of the absolute 15-25 frequency Lo ge /2-K, 
if we restrict ourselves to the experimental range of pow- 
ers (0.1 W < P < 0.5 W). This does not contradict the 
uncertainty of the line shape model of 20 Hz given in Q . 
Moreover, in the derivation of the possible small tempo- 
ral drift of the 15-25 frequency Q this error cancels out. 
Still, for a further improvement of the 15-25 frequency 
measurement it is highly desirable to reduce the laser in- 
tensity and thus reduce the error and the uncertainty of 
the extrapolation. 

Secondly, the calculations show, that due to the pho- 
toionization, the velocity distribution of atoms escaping 



the nozzle significantly differs from the distribution of 
atoms contributing to the signal. For the correction of 
the second order Doppler effect it is more physical to 
model the initial velocity distribution in the nozzle by a 
Maxwellian distribution (see Eq. Q), which is dynam- 
ically modified by the ionization process of atoms with 
long interaction times, instead of modelling an effective 
velocity distribution at the detector. On the other hand, 
the magnitude of the Doppler effect, which has been one 
of the biggest concerns in the last measurements 0, 
is already on the level of —20 Hz for a delay time of 
At = 1210 /is. This opens an opportunity to use a sim- 
plified line shape model to correct for this effect. 

We have performed an evaluation of the whole data 
set obtained in 2003 measurement [j| using such a sim- 
plified approach. The 15-25 spectral lines detected with 
the delays At > 1210 fis have been fitted by Lorentzian 
functions. The frequencies of the line centers have been 
corrected for the second-order Doppler effect by adding 
the calculated frequency shift. The result of this analysis 
was consistent with the result of the full line shape model 
9] within the uncertainty of the model. The shift should 
be sensitive to the initial velocity distribution which can 
differ from due to e.g. geometrical factors, unperfect 
thermalization or collisions. The delayed detection mod- 
ifying the velocity distribution (J5J drastically reduces the 
mean Doppler shift as well as its sensitivity to the initial 
distribution. We tested numerically, that the substitu- 
tion of the v 3 dependency by v A results in a mean 
Doppler shift of —23 Hz for At = 1210 /its, which is only 
3 Hz aside the previous value of —20 Hz. The difference 
decreases for increasing At and the corresponding un- 
certainty in the determination the 15-25 frequency can 
be reduced to a few hertz. Thus, using such a simplified 
model we can correct the second-order order Doppler ef- 
fect with an uncertainty on the level of a few parts in 
10 15 . 

The last point concerns the influence of the freezing 
nozzle on the observed absorption line centers, and re- 
sulting effects on the absolute frequency, which is ob- 
tained by the extrapolation of these line centers to zero 
laser power. Since the nozzle temperature is lower than 
its melting point of 14 K, molecular hydrogen effectively 
freezes on the nozzle walls and reduces its diameter. In 
the experiment, we had to melt molecular hydrogen ap- 
proximately once per 20 min to avoid significant losses 
of laser power and atomic hydrogen flux in the 243 nm 
cavity. Our simulation allows for different spatial seed- 
ings of starting points of the atomic trajectories at the 
diaphragms, which can be used to model a reduced noz- 
zle diameter or a misalignment of the cavity mode with 
respect to the atomic beam axis. 

We have performed a simulation of a whole day of mea- 
surement recording the 1210 fis delayed spectra at 100 
different intensities, while varying the nozzle radius ran- 
domly in the range from rui down to wq. The resulting 
line center fits as a function of laser power per direction 
are shown in Fig. |3p. Figure [S^ results from the same 
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FIG. 5: Monte-Carlo simulation: (a) - the frequency of the 
15- 25 transition vs. laser power for an ensemble of 10 000 
atoms for the completely opened nozzle. The residual scatter 
of the data can be explained by the stochastic nature of the 
Monte-Carlo simulation, (b) - the same simulation, but the 
nozzle radius is randomly varied in the range {wo,rDi}. We 
fit the data at low intensity limit by linear functions (dashed 
tilted lines) for visualization of non-linearity of the frequency 
shift. 



simulation but with a permanently free nozzle. An un- 
controlled freezing leads to an excess scatter of the line 
centers, increasing for larger exciting powers. Specifi- 
cally, a smaller radius at the entrance diaphragm results 
in an upwards shift in line center frequency, the shift in- 
creasing with excitation power. This is consistent with 
the observation in Fig. f4J:, where a small nozzle radius 
results in a larger value for fc s hift- As a result, a linear fit 
in the range of 0.1-0.5 W and an extrapolation to zero 
intensity gives different values for the absolute frequency 
depending on the amount of nozzle freezing. Typically 
the scatter is on the level of a few tens of Hertz, which 
can only partly explain the scatter observed in the exper- 
iments [fj, 13 • Other contributions can originate in un- 
certainties of the power axis and misalignment of the 243 
nm enhancement cavity. In future measurements we will 
aim at increasing the nozzle diameter to avoid restriction 
of the atomic trajectories, working at lower excitation 
powers, and improving the stability of the enhancement 
cavity including power stabilization. 

The beam nature of the 15-25 spectroscopy of atomic 
hydrogen imposes certain restrictions on the experimen- 



tal conditions like high intensity of the excitation field, 
impossibility of further cooling of the beam, short inter- 
action time and freezing of the nozzle. The enumerated 
effects result in a number of systematic shifts and cor- 
responding loss of accuracy. Note, that all optical spec- 
troscopic experiments demonstrating the accuracy better 
than 10 -14 deal with cold ions or atoms in traps (see e.g. 
EH El)- All the systematic shifts of the 15-25 transi- 
tion discussed in the paper would be significantly reduced 
or eliminated if the hydrogen atoms were slowed down by 
laser cooling or if one used hydrogen-like ions that could 
be held in a trap E3 ■ 

Based on considerations of this Section we can con- 
clude, that perturbation of the mode of the enhancement 
cavity results in a frequency shift via the intensity related 
shift and broadenings. Though the cavity is not confo- 
cal which ensures the suppression factor over 500 for the 
TEMoi mode, the intra-cavity diaphragms can influence 
the structure of the TEMnn cavity mode. Monte-Carlo 
simulations performed for impure TEMnn mode shows 
that the results presented in the article remain valid if 
the power falling into the TEMni mode does not exceed 
5%. Thus, special attention should be paid for inten- 
sity stabilization, coupling and alignment of the cavity 
to avoid any additional non resonant light. After taking 
the issues discussed in the paper into account it seems 
realistic to reduce the uncertainty of the determination 
of the 15-25 frequency to the level of a few parts in 
10 15 . 



V. CONCLUSIONS 

In this article we discussed the transient excitation dy- 
namics of the 15-25 transition of atomic hydrogen un- 
der irradiation of two counter-propagating 243 nm light 
beams. The main focus of the present analysis has 
been the inclusion of photoionization of the 25 state, 
power broadening and the inhomogeneous AC Stark in 
the quantum dynamics and the detailed study of their 
influence on the transition frequency and its line width. 

Though photoionization would only insignificantly 
change the line width of each individual atom in the 
experimental range of excitation powers and interaction 
times in a laser field of constant intensity, the preferred 
ionization of slow 25 atoms results in a substantial broad- 
ening of the observed collective spectrum. The mod- 
elling shows, that the resulting intensity-dependent fre- 
quency shift is mainly caused by the dynamic Stark 
shift, while the intensity-dependent line broadening re- 
sults from power broadening, the dynamic Stark shift in 
a spatially inhomogeneous laser field and the influence of 
photoionization. 

For our simulation we used the values of the transition 
matrix elements, the dynamic Stark shift coefficients and 
the ionization rates re-derived in The results of the 
simulation show a fair agreement with the experimentally 
observed ones. We can confirm the value of the difference 
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of the dynamic Stark shift coefficients /3 ac (25) — /3 ac (15) 
with an accuracy of 15% using the analysis of the 15-25 
frequency shift. The analysis of the line broadening also 
shows a good agreement with the experimental values 
within 18%. Due to a complex structure of the intensity- 
dependent line broadening it does not allow to make def- 
inite conclusions about the ionization coefficient, though 
indicating that our treatment of the excitation process in 
the thermal beam is correct. 

The uncertainty of the determination of the 15-25 
frequency is determined not only by the experimental 
accuracy, but also by the uncertainty of the line shape 
model and the correction of the dynamic Stark shift. 
Taking photoionization into account allows to build a 
more accurate line shape model. Factors such as non- 
linearity of the line shift and geometrical factors con- 
tribute on a level comparable in magnitude with the 



second order Doppler effect, which has been of major 
concern in previous studies 0,01 • The Monte-Carlo ap- 
proach presented in this paper allows to perform accurate 
quantitative evaluations of all these effects and make an 
accurate evaluation of the final uncertainty. 
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